Chaotic dynamics of a glaciohydraulic model

A model subglacial drainage system, coupled to an ice-dammed reservoir that receives a time-varying meltwater input, is described and analysed. In numerical experiments an ice-marginal lake drains through a subglacial channel, producing periodic floods, and fills with meltwater at a rate dependent on air temperature, which varies seasonally with a peak value of Tm. The analysis reveals regions of Tm parameter space corresponding to ‘mode locking’, where flood repeat time is independent of Tm; resonance, where decreasing Tm counter-intuitively increases flood size; and chaotic dynamics, where flood cycles are sensitive to initial conditions, never repeat and exhibit phase-space mixing. Bifurcations associated with abrupt changes in flood size and timing within the year separate these regions. This is the first time these complex dynamics have been displayed by a glaciohydraulic model and these findings have implications for our understanding of ice-marginal lakes, moulins and subglacial lakes

Antarctic surface hydrology and impacts on ice-sheet mass balance

Melting is pervasive along the ice surrounding Antarctica. On the surface of the grounded ice sheet and floating ice shelves, extensive networks of lakes, streams and rivers both store and transport water. As melting increases with a warming climate, the surface hydrology of Antarctica in some regions could resemble Greenland’s present-day ablation and percolation zones. Drawing on observations of widespread Antarctica surface water and decades of study in Greenland, we consider three modes by which meltwater could impact Antarctic mass balance: increased runoff, meltwater injection to the bed, and meltwater-induced ice-shelf fracture, all of which may contribute to future ice sheet mass loss from Antarctica

Modelling channelized surface drainage of supraglacial lakes

Supraglacial lakes can drain to the bed of ice sheets, affecting ice dynamics, or over their surface, relocating surface water. Focusing on surface drainage, we first discuss observations of lake drainage. In particular, for the first time, lakes are observed to drain >70 km across the Nivlisen ice shelf, East Antarctica. Inspired by these observations, we develop a model of lake drainage through a channel that incises into an ice-sheet surface by frictional heat dissipated in the flow. Modelled lake drainage can be stable or unstable. During stable drainage, the rate of lake-level drawdown exceeds the rate of channel incision, so discharge from the lake decreases with time; this can prevent the lake from emptying completely. During unstable drainage, discharge grows unstably with time and always empties the lake. Model lakes are more prone to drain unstably when the initial lake area, the lake input and the channel slope are larger. These parameters will vary during atmospheric-warming-induced ablation-area expansion, hence the mechanisms revealed by our analysis can influence the dynamic response of ice sheets to warming through their impact on surface-water routing and storage

Geothermal heat flux is the dominant source of uncertainty in englacial-temperature-based dating of ice rise formation

Ice rises are areas of locally grounded, slow-moving ice adjacent to floating ice shelves. Temperature profiles measured through ice rises contain information regarding changes to their dynamic evolution and external forcings, such as past surface temperatures, past accumulation rates and geothermal heat flux. While previous work has used borehole temperature–depth measurements to infer one or two such parameters, there has been no systematic investigation of parameter sensitivity to the interplay of multiple external forcings and dynamic changes. A one-dimensional vertical heat flow forward model developed here examines how changing forcings affect temperature profiles. Further, using both synthetic data and previous measurements from the Crary Ice Rise in Antarctica, we use our model in a Markov chain Monte Carlo inversion to demonstrate that this method has potential as a useful dating technique that can be implemented at ice rises across Antarctica. However, we also highlight the non-uniqueness of previous ice rise formation dating based on temperature profiles, showing that using nominal values for forcing parameters, without taking into account their realistic uncertainties, can lead to underestimation of dating uncertainty. In particular, geothermal heat flux represents the dominant source of uncertainty in ice rise age estimation. For instance, in Crary Ice Rise higher heat flux values (i.e. about 90 mW m−2) yield grounding timing of 1400 ± 800 years, whereas lower heat flux of around 60 mW m−2 implies earlier ice rise formation and lower uncertainties in the ice rise age estimations (500 ± 250 years). We discuss the utility of this method in choosing future ice drilling sites and conclude that integrating this technique with other indirect dating methods can provide useful constraints on past forcings and changing boundary conditions from in situ temperature–depth measurements.</p

Far-field optical microscope with nanometer-scale resolution based on in-plane surface plasmon imaging

A new far-field optical microscopy technique capable of reaching nanometer-scale resolution has been developed recently using the in-plane image magnification by surface plasmon polaritons. This microscopy is based on the optical properties of a metal-dielectric interface that may, in principle, provide extremely large values of the effective refractive index n up to 100-1000 as seen by the surface plasmons. Thus, the theoretical diffraction limit on resolution becomes lambda/2n, and falls into the nanometer-scale range. The experimental realization of the microscope has demonstrated the optical resolution better than 50 nm for 502 nm illumination wavelength. However, the theory of such surface plasmon-based far-field microscope presented so far gives an oversimplified picture of its operation. For example, the imaginary part of the metal dielectric constant severely limits the surface-plasmon propagation and the shortest attainable wavelength in most cases, which in turn limits the microscope magnification. Here I describe how this limitation has been overcome in the experiment, and analyze the practical limits on the surface plasmon microscope resolution. In addition, I present more experimental results, which strongly support the conclusion of extremely high spatial resolution of the surface plasmon microscope.Comment: 23 pages, 9 figures, will be published in the topical issue on Nanostructured Optical Metamaterials of the Journal of Optics A: Pure and Applied Optics, Manuscript revised in response to referees comment

Late Holocene ice-flow reconfiguration in the Weddell Sea sector of West Antarctica

Here we report Late Holocene ice sheet and grounding-line changes to the Weddell Sea sector of West Antarctica. Internal radio-echo layering within the Bungenstock Ice Rise, which comprises very slow-flowing ice separating the fast-flowing Institute and Möller ice streams, reveals ice deformed by former enhanced flow, overlain by un-deformed ice. The ice-rise surface is traversed by surface lineations explicable as diffuse ice-flow generated stripes, which thus capture the direction of flow immediately prior to the creation of the ice rise. The arrangement of internal layers can be explained by adjustment to the flow path of the Institute Ice Stream, during either a phase of ice sheet retreat not longer than ∼4000 years ago or by wholesale expansion of the grounding-line from an already retreated situation not sooner than ∼400 years ago. Some combination of these events, involving uplift of the ice rise bed during ice stream retreat and reorganisation, is also possible. Whichever the case, the implication is that the ice sheet upstream of the Bungenstock Ice Rise, which currently grounds over a >1.5 km deep basin has been, and therefore may be, susceptible to significant change

Major ice‐sheet change in the Weddell Sector of West Antarctica over the last 5000 years

Until recently, little was known about the Weddell Sea sector of the West Antarctic Ice Sheet. In the last 10 years, a variety of expeditions and numerical modelling experiments have improved knowledge of its glaciology, glacial geology, and tectonic setting. Two of the sector's largest ice streams rest on a steep reverse‐sloping bed yet, despite being vulnerable to change, satellite observations show contemporary stability. There is clear evidence for major ice‐sheet reconfiguration in the last few thousand years, however. Knowing precisely how the ice sheet has changed in the past, and when, would allow us to better understand whether it is now at risk. Two competing hypotheses have been established for this glacial history. In one, the ice sheet retreated and thinned progressively from its Last Glacial Maximum position. Retreat stopped at, or very near, the present position in the Late Holocene. Alternatively, in the Late Holocene the ice sheet retreated significantly upstream of the present grounding line, and then advanced to the present location due to glacial isostatic adjustment, and ice‐shelf and ice rise buttressing. Both hypotheses point to data and theory in their support, yet neither has been unequivocally tested or falsified. Here, we review geophysical evidence to determine how each hypothesis has been formed, where there are inconsistencies in the respective glacial histories, how they may be tested or reconciled, and what new evidence is required to reach a common model for the Late Holocene ice sheet history of the Weddell Sea sector of West Antarctica

Constraining recent ice flow history at Korff Ice Rise, West Antarctica, using radar and seismic measurements of ice fabric

The crystal orientation fabric of ice reflects its flow history, information which is required to better constrain projections of future ice sheet behavior. Here we present a novel combination of polarimetric phase‐sensitive radar and seismic anisotropy measurements to provide independent and consistent constraints on ice fabric at Korff Ice Rise, within the Weddell Sea sector of West Antarctica. The nature and depth distribution of fabric in the ice column is constrained using the azimuthal variation in (1) the received power anomaly and phase difference of polarimetric vertical radar soundings and (2) seismic velocities and shear wave splitting measurements. Radar and seismic observations are modeled separately to determine the nature and strength of fabric within the ice column. Both methods indicate ice fabric above 200‐m depth which is consistent with present‐day ice‐divide flow. However, both measurements also indicate an oblique girdle fabric below 230‐m depth within the ice column, inconsistent with steady state divide flow. Our interpretation is that this deeper fabric is a remnant fabric from a previous episode of flow, which is currently being overwritten by ongoing fabric development associated with the present‐day flow regime. The preexisting fabric is consistent with ice flow from the south prior to ice‐divide formation, in agreement with models of Holocene ice sheet evolution. These findings apply new constraints to the flow history at Korff Ice Rise prior to divide formation and demonstrate the capacity of radar and seismic measurements to map fabric and thus constrain past ice flow